Temperature Dependence and Screening Models in Quantum Cascade Structures 
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Different screening models in quantum cascade lasers are compared by calculating the contribution 
of intra-subband impurity scattering to the optical linewidth as a function of temperature. We find 
a strong impact of impurity scattering which is increasing substantially with temperature. A simple 
isotropic bulk screening model works well if the screening length is of the order or longer than the 
period length of the cascade structure. 



I. INTRODUCTION 

The performance of terahertz (THz) quantum cascade 
lasers 1 (QCLs) has improved fast, now working up to 186 
K— , and with frequencies down to 1.2 THz 3 without the 
use of magnetic fields. A main issue of current research 
is to reach higher operation temperatures^. This requires 
a good understanding of the temperature dependence 
for all components contributing to laser operation. Of- 
ten, phonon interaction is addressed when temperature- 
dependent effects in lasers are discussed^. However, as 
previously suggested^, there is also a strong impact on 
operation by screening, which strongly depends of the 
electron temperature. 

The fast temperature degradation of THz QCLs has 
previously been investigated through simulation. Indjin 
et.al^ found a strongly increasing scattering rate with 
temperature from the upper to the lower laser state by 
optical phonon emission, which lowers the population 
inversion. This effect has also been quantified by Ji- 
rauscheck and Lugli 9 and Cao et.al. 6 . One approach to 
remove this effect is by applying a strong static mag- 
netic fielc£ and also, increasing the optical phonon en- 
ergy by changing to other material systems such as gal- 
lium nitride is expected to lower this effect. However, 
this thermally activated phonon emission cannot explain 
the empirical relation that the maximum temperature is 
approximately limited by the photon energy 4 
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since if this was the dominant degradation effect, long- 
wavelength lasers should have slightly better temperature 
performance. 

A second phonon effect discussed in this context is 
thermal backfilling, i.e., electrons absorb optical phonons 
and are reinjected to the lower laser state. This process 
has been theoretically shown to play a minor role^ and, 
also, the double phonon resonance did not improve the 
temperature performance of the device^. 

A simplified relation for the gain at resonance is [see 
e.g. Eq. (1) in Ref . [n| 
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where A/ is the population inversion and r spect is the 



FIG. 1: Illustration of two, partly compensating, tempera- 
ture dependent broadening effects, (a) As the temperature is 
increased, the electrons populate states with higher momen- 
tum, (b) Due to the constant density of states, the number 
of final states are constant. However, regarding elastic intra- 
band scattering, for large momentum, the average momentum 
transfer is larger than for small momentum, (c) The scat- 
tering matrix element as a function of momentum transfer. 
Although increasing with temperature, it is decreasing with 
momentum transfer. 



spectral width of the transition due to scattering pro- 
cesses. Lasing operation requires <?i as to overcome the 
losses, and thus Eq. (JTJ) is immediately satisfied if ei- 
ther A/ cx X/T, or r spcc t cx T. (Of course a combi- 
nation between these together with the frequency de- 
pendence of losses is more realistic for a full quanti- 
tative understanding). Most previous studies have fo- 
cused on A/, while r spcct is also expected to increase 
with temperature^. Similarly, Li et.al^- found recently 
that even though thermally activated phonon emission 
strongly reduced the population inversion with tempera- 
ture, life-time broadening had to be included in the model 
in order to correctly estimate the experimental maximum 
lasing temperature. In this manuscript we investigate the 
temperature dependence of impurity scattering at ion- 
ized dopants, which are present in all lasers considered. 
This scattering process is strongly affected by screening, 
which, as a many-body effect, is difficult to treat numer- 
ically and approximations are necessary. In this work we 
investigate different approximations regarding screening 
and their impact on the scattering strength. 

The impact of temperature on this impurity scatter- 
ing process is dominated by two competing processes: 
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On the one hand, the scattering strength is increasing 
with temperature due to the decrease of screening. There 
is, however, a compensating effect: At low temperatures 
electrons mainly occupy the lower parts of the subband, 
where other momentum states at the same energy only 
differ by small momentum. Small momentum transfer 
give rise to a large impurity scattering matrix element 
and therefore strong scattering. At higher temperatures, 
electrons populate the upper part of the subbands, where 
the average momentum transfer is larger, giving rise to a 
smaller scattering rate, see Fig. [1] The case of screening 
in electron-electron scattering is even more difficult due 
to the necessary scattering partner. 

The manuscript is organized as follows: First, we 
present the theory used to treat screening based on the 
random phase approximation. Then we show results for 
the scattering strength including also screening from an 
infinite periodic structure. After these findings we dis- 
cuss the relation between the scattering strength and the 
actual spectral linewidth of the transition. 



II. THEORY 

The method used here is based on the work of Lee and 
Gailbraithi^ regarding the screening and Andoi^ con- 
cerning impurity scattering and broadening. 

In the multi-subband structures studied here, the 
Dyson-like equation governing the screened Coulomb 
scattering matrix elements, Wijki(q), is of the form, 

Wijki(q) = Vijki(q) + ^2 v mm(q)^mn.(q)W mn ki(q) (3) 

mn 

where Vijki{q) is the unscreened, or bare, Coulomb ma- 
trix element, n m „(g) is the polarization function, q is 
the magnitude of the in-plane momentum vector, and 
the indices refer to different subbands. The unscreened 
Coulomb matrix element is 



Vi jk i (q) 
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where e < is the electron charge, A the in-plane area of 
the device, eo is the vacuum permittivity, and s r is the 
relative permittivity of the background material. The 
form factor is 

F ijkl {q) = JdzJ dz'i,*(z)^(z)e-«\*-*yUz')Mz')- 

(5) 

Solving these double integrals can be strongly simplified 
by the use of the scheme presented in the latter half of 
Sec. 2B in Ref. GJ. 

The important quantity in this context is the polar- 
ization function, where most approximations are usually 
done. In this work we will use the (static) random-phase 
approximation (RPA), 
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where f m ,k (£ m ,k) is the occupation (energy) of state k 
in subband m, and the factor 2 is for spin. Throughout 
this work we approximate the distributions / with a ther- 
malized Fermi-Dirac distribution using the total electron 
densities densities ni per subband obtained from our de- 
vice simulation under operating conditions. If not stated 
otherwise we use temperature-independent n, (evaluated 
at I00K) and vary all band temperatures equally in order 
to focus on the screening effect. 

The expression, similar to Eq. ([3]), for determining the 
screened impurity scattering matrix element is 
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where ^™ p (g) is the unscreened impurity scattering ma- 
trix element, 
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where z- lmp is the spatial location of the impurity in the 
growth direction. 

The linewidth of the transition at wave vector k, is es- 
timated by the average of the intra-band scattering rates 
of the two lasing states, 

T(k) = nJ2 <Kl(3) + Vu(q)) s . c . <^k +q - E k ) , (9) 



where u and / are the indices of the subbands between 
which the lasing transition occurs and (■■■) sc is the 
average over different scattering configurations. Inter- 
subband scattering is weaker and neglected here. How- 
ever, the temperature dependence of screening is ex- 
pected to follow the same trend in both types of scat- 
tering processes. The average scattering strength related 
to the lasing transition is then given by 



£ k r(k)|/ U|k -/ t , k | 

J2k \ fu,k ~ fl,k\ 
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In simple models, this scattering strength can be directly 
related to the spectral linewidth of the transition. How- 
ever, this has to be taken with care as discussed in Section 

m 

In this work, we compare three different screening 
models: (i) The full RPA model [Eqs. 1(515)1]. where we 
include all states from N periods of the QCL structure 
surrounding the upper and lower laser states. This is 
denominated by RPA-A^. (ii) The long wavelength limit 
(q — > 0), where 
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Ei — Ej 



ifi^i, (11) 



with rij = 2/AJ2k /i,k=ch as used, e.g., in the work of Lii 
and Cao^, and (iii) the simple isotropic screening model 
where the screened impurity matrix element is obtained 
by the replacement q — > \J q 2 + A 2 in the right-hand side 
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FIG. 2: Scattering strength due to intra-band impurity 
scattering for the three different screening models and the 
isotropic Debye approximation. For all screening models, the 
scattering is increasing with temperature due to less screen- 
ing. Both the isotropic screening model and the long wave- 
length limit model slightly under-estimate the scattering (by 
over estimating the screening) compared to the RPA result, 
but are, however, excellent approximations in this situation. 



of Eq. ([5]) , used by us earlier—. A is the inverse screening 
length, 

A 2 = ^^£ (12) 

which, in the Debye approximation (valid if fc^T > /_*), 
becomes 

A ° cb ^ " e r e Q k B r [16 > 

where n^D is the average electron density in the struc- 
ture and \x is the chemical potential of the corresponding 
three-dimensional electron gas. 



III. RESULTS 
A. Scattering Strength 

We focus on the THz QCL presented in Ref. [Tt] using 
identical parameters as in Ref. [D. Five states per period 
and three periods are included in the simulation. 

The scattering strength due to impurity scattering is 
depicted in Fig. [21 The scattering, in all approximations, 
is increasing with temperature due to the reduction of 
screening with temperature. The reduction of scatter- 
ing matrix elements with momentum does not depend as 
strongly on temperature, as transitions with small mo- 
mentum transfer are always present. 

In Fig. [3] the scattering strength is depicted for vary- 
ing subband populations according to the results from 



FIG. 3: Scattering strength calculated by RPA-3 for different 
choices of band temperature and densities. The temperature 
dependence of subband populations m (T) hardly changes the 
results, while increasing the temperature of the lasing sub- 
bands by 100 K lowers the scattering strength. 

the non-equilibrium Green's function transport simula- 
tion presented in Ref. 0- The varying subband popula- 
tions have a negligible effect on the scattering strength, 
even though the population in, e.g., the lower lasing sub- 
band increases by a factor of 2 over the studied temper- 
ature interval. 

The electronic temperature in the heavily populated 
injector subband has been measured to be only slightly 
higher than the lattice temperature, while the less pop- 
ulated lasing subbands have temperatures of the order 
100 K higher than the lattice^. It is therefore of im- 
portance to calculate the scattering strength in this non- 
equilibrium situation. Fig. [3] shows that heating of the 
laser level populations decreases the scattering. However, 
the overall trend of increasing scattering with tempera- 
ture persists. As demonstrated in Fig. 0] the screened 
interaction T(k) is more affected by the temperature of 
the injector level than that of the laser level. This shows 
that screening is dominated by the majority of electrons 
in the injector level and accordingly their temperature is 
the most relevant. The reduction of scattering strength 
by heating of the laser levels, observed in Fig. [3j is at- 
tributed to the q dependence of the scattering rate, where 
higher values of q are relevant for heated distributions. 

B. Dependence on the number of periods 

For the data shown so far, only a finite set of periods 
are included in the screening models, which is a com- 
mon approximation 15,16 . For comparison, we study now 
the infinitely periodic system, reflecting the large num- 
ber of periods in real QCL structures. First, we need to 
make an approximation. Up until now we have taken all 
combinations of states into account when studying the 
electron-electron matrix elements, Vijki- We will now ne- 
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FIG. 4: T(k) (full lines) and k\f u ,k — fi,k\ (dashed lines, where 
the prefactor k mimics the density of states in the subband). 
In the first case all subband temperatures are set to 50 K. In 
the second case the lasing subbands are 100 K warmer, and 
hence, the individual rates F(k) are larger due to a slightly 
reduced screening. However, the electrons in lasing subbands 
affect the weight function fc|/u,k — /j,k| which is shifted to 
higher momentum states and therefore give a smaller average 
scattering strength. In contrast, if the injector subband is 
100 K warmer (third case) there is even less screening and 
stronger scattering. This shows that the injector subband 
(with 59 % of the population) is more important for screening 
than the lasing subbands (with 30 % of the population). If all 
subbands are set to 150 K (fourth case) scattering becomes 
even stronger so that the average scattering strength increases 
taking into account the shift to higher momenta. 

gleet all matrix elements where i and j do not belong 
to the same period and/or k and I do not belong to 
the same period. This would correspond to restricting 
to only intra-period scattering while keeping both intra- 
and inter-period interaction. Taking three periods into 
account, this approximation has only a minor effect on 
the scattering strength, see Fig. [5] The notation can then 
be changed to Viju-.h, where the first four indices denote 
subbands within a single period and h is an integer de- 
noting the distance in periods between the ?j-pair and 
the fcZ-pair. h — corresponds to that all subbands are 
in the same period, h = 1 to the /cZ-pair being in the first 
period to the right of ij, and so on. 

The Dyson equation, Eq. ([3]), in this new notation, 
becomes 

WijkLh = VijkLh + ^ Vijnm,h'^-mnW mn kl,h-h' (14) 
ran.h' 

omitting the wave vector dependence. This convolution 
in h decouples in Fourier space. Therefore, we introduce 
the transformation, 

': r (15) 

2tt J -tv 



FIG. 5: Scattering strength for different screening models. As 
an estimate of the error introduced when restricting to intra- 
period scattering RPA-3 with and without inter-period scat- 
tering is plotted. At low temperatures, screening is efficient 
and all RPA models give the same results. At higher temper- 
atures, the periodic model approaches the isotropic model. 



where Eq. (fT4"|) becomes 

W ijkllK = V m , K + £ V ijnm>K U mn W mnkl<K . (16) 
mn 

For large h, the wavefunctions of the two pairs do not 
overlap spatially [see Eq. ([5])] and the matrix elements 
become particularly simple. If H is the largest h where 
the wavefunctions overlap, the matrix elements becomes 
V ijklih (q) = e~^ h - H ^ d V l]kl ,H(q) for h > H where d is 
the length of one laser period. This exponential h de- 
pendence makes the sum over h in the Fourier transform, 
Eq. (fT3|) . simple and can be easily be evaluated to infin- 
ity via a geometric sum. In computations, the symmetry 
Vijki-h = Vkuj--h will be used. Eq. (fT6|) can be solved for 
a finite number of n, and the screened impurity matrix 
element can then be calculated for infinite system. 

The results can be seen in Fig. [5] First, we notice 
that including only one period in screening strongly over- 
estimates scattering. Also, when including three periods, 
restricting to only intra-period scattering underestimates 
screening slightly. Finally, we see that electrons in pe- 
riods further away than in the nearest neighboring peri- 
ods, i.e., the difference between RPA-oo and RPA-3 with 
intra-period scattering, give a moderate effect on screen- 
ing at higher temperatures. 

This periodic formulation does not only give the possi- 
bility to investigate the screening contribution of electron 
in arbitrary periods, but also, the scattering impact on 
the lasing transition from impurities in any period can 
easily be calculated, see Fig. [6] We see that scattering 
off dopants from the same periods is the strongest one, 
but also, at higher temperatures, dopants in neighboring 
periods give a substantial contribution to the scattering. 
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T= 10 K — — 




FIG. 6: Scattering strength due to impurity layers in neigh- 
boring periods. At low temperatures screening is efficient and 
scattering is spatially local. When temperature is increased, 
scattering at the impurities in the same period is increased, 
but also, scattering at impurities in neighboring periods be- 
comes important. 



C. Mid-Infrared Lasers 

Up until now, only the THz device presented in Ref. [ij] 
has been investigated. QCLs are divided into two differ- 
ent types: THz lasers which operates at frequencies be- 
low the optical phonon resonance, and mid-infrared (mid- 
IR) laser operating at frequencies above. Typically, mid- 
IR lasers are more heavily doped and, in order to reach 
higher optical frequencies, have more and narrower quan- 
tum wells. The period in these lasers can often clearly 
be divided into two spatial regions: one active region, to 
which the lasing subbands are confined, often only three 
quantum wells, and an injector region where the dopants 
are located. This spatial separation between the dopants 
and lasing subbands together with the higher electron 
concentration contribute to effective screening that low- 
ers the scattering impact from impurities on the lasing 
transition. Also, other scattering mechanisms, such as 
interface roughness become much more important due to 
the larger conduction band offset. 

The results from our screening model can be seen in 
Fig. [7] and [5J where the laser presented in Ref. [l^ has 
been used. First, the impact of impurity scattering on the 
transition is much smaller for reasons mentioned above, 
but also, the simple isotropic screening model gives a 
larger error. This can be understood by comparing the 
screening lengths. The two lasers in this study have sim- 
ilar period, but the mid-IR device is much more doped 
and therefore has a shorter screening length, see Tab. fl] 
If the screening length is of the order of the laser period 
or longer, all electrons in the structure will contribute 
to screening and the isotropic screening model will be a 
good approximation. Then, the spatial location of the 
electrons does not matter since all electrons contribute 



TABLE I: Comparison between typical THz and mid-IR de- 
vices with index h. Note that the high doping in the mid- 
IR device gives a degenerate electron gas at low tempera- 
tures resulting in an almost temperature independent screen- 
ing length. In order to get a screening length in the mid-IR 
device of half the period length, a temperature of 720 K is 
required. For the THz device, the electron gas is well approx- 
imated by a non-degenerate gas resulting in Debye screening, 
see Eq. (JT^I 

THz Device 17 Mid-IR Device 19 

d 55.4 nm 45.3 nm 

n ZD 4.06 x 10 15 cm -3 8.6 x 10 16 cm -3 

1/A(50K) 27.9 nm 7.46 nm 

1/A(200 K) 55.0 nm 10.2 nm 
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50 100 150 200 250 300 
T[K) 

FIG. 7: Scattering strength for different approximations for 
the mid-IR device of Ref. Il9l . The curves for the RPA with 
2 and infinite periods (as well as the approximation q — » 0) 
fall together, showing that the strong screening in this heavily 
doped structure restricts the interaction to the two neighbor- 
ing periods. Changing the subband population according to 
a simulation based on non-equilibrium Green's functions af- 
fects the screening more in this laser compared to the THz 
structure due to the strong, local screening. At temperatures 
above 300 K the different scattering strengths increase almost 
linearly. Here, the RPA-oo and the isotropic screening model 
exhibit an almost constant difference, thus the relative error 
becomes smaller. The RPA-2 starts to deviate from RPA-oo 
around 700 K, where 1/A = 22.5 nm. 



to screening. In the mid-IR device, the screening length 
is much shorter than the period and the isotropic model 
screening fails. Also, the screening is more affected by 
temperature dependent subband populations, since the 
screening is strongly depending on the electron concen- 
tration in the direct vicinity of the dopants. Thus the 
isotropic screening model is an excellent approximation 
in the lightly doped THz quantum cascade lasers but is 
questionable for highly doped mid-IR devices, where the 
screening length is of the order of the individual layers. 
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FIG. 8: Scattering strength due to impurity layers in neigh- 
boring periods for the mid-IR device of Ref.—. As the ac- 
tive region where the lasing transition occurs is sandwiched 
between two injector regions where the dopants are located, 
the scattering contributions from each side is approximately 
equal. Due to the heavy doping, impurity layers in next neigh- 
boring periods are completely screened, also at high temper- 
atures. 



IV. RELATION BETWEEN SCATTERING 
STRENGTH AND SPECTRAL LINEWIDTH 

Calculating the spectral width of the gain transition 
in QCLs is quite non-trivial due to correlation effects^. 
A simple closed expression for the linewidth including 
correlation effects can be obtained by simply replacing 
Eq. © by 

r corr (k) =irJ2 (\Vuu(q) ~ Va(q)\ 2 ) s c <5(£ k+q - E k ) , 
q 

(17) 

see Ref. [jj]. This form clearly shows the effect of corre- 
lation effects, namely that if the scattering environment 
of the two lasing states are similar, the optical linewidth 
is reduced below the sum of the lifetime induced widths 
of the two respective states. Also, the estimated width 
of the transition is now a measure of the difference in 
scattering strength and therefore does no longer directly 
reflect the scattering strength itself. The result is shown 
in Fig.[5]for different approximations. Again the isotropic 
screening model allows for a quantitative description in 
comparison with the more advanced RPA-approaches. 
Compared to Fig. [2l the estimated width is strongly re- 
duced and decreasing with temperature except for the 
low temperature regime. This difference can be under- 
stood by the fact that the impurity scattering matrix 
elements for q — ► 0, which are most affected by screen- 
ing, widely cancel each other in Eq. (jTTJ) , as outlined in 
the Appendix. 

However, one has to bear in mind, that Eq. (jTTJ) seems 
to under-estimate the optical linewidth compared to a 
full calculation, see Ref. 11. This goes in hand with the 



FIG. 9: Scattering strength due to intra-band impurity scat- 
tering inclu ding correlation effects by Eq. (|17[) f° r the THz 
laser of Ref. Ha 



findings of Ref. 7, where a full calculation based on the 
isotropic, bulk screening model found an increase of the 
Full Width at Half Maximum from 2.7 meV at 10 K to 
5.5 meV at 80 K, which is between the results of the total 
scattering strength from Fig. [3] and the correlated scat- 
tering strength from Fig. [HI Thus, both the expressions 
© and (JTTJ) should be taken with care for the determi- 
nation of the width of the gain peak. 

Finally, the averaging in Eq. (fTU|) is over-estimating the 
linewidth since the averaged gain peak has a linewidth 
smaller than the average linewidth. However, this effect 
is expected to be weak compared to the correlation effects 
addressed aboveJ^ 



V. CONCLUSION 

In conclusion, the decrease of screening causes an in- 
crease of scattering strength in QCLs with temperature. 
If the screening length is of the order of the laser pe- 
riod or longer, the isotropic, bulk screening model is an 
excellent approximation. In this situation, which is com- 
mon for typical THz QCLs, the distribution of electrons 
to different subbands is much less important than the 
average subband temperature. For high-doped mid-IR 
laser structures, where the screening length is typically 
of the order of the layers, a microscopic RPA model is 
indispensable for a reliable calculation of screening. The 
temperature dependence of the linewidth is difficult to es- 
timate due to correlation effect, and a full self-consistent 
transport calculation is needed to resolve this matter. 

The authors thank C. Weber, C.-O. Almbladh and 
M. P. von Friesen for helpful discussions and gratefully 
acknowledge financial support from the Swedish Research 
Council (VR). 
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APPENDIX A: IMPACT OF SCREENING ON 
CORRELATIONS IN SCATTERING MATRIX 
ELEMENTS 

Screening affects the interaction over larger distances, 
corresponding to small wave vectors. In the limit q — *• 0, 
the diagonal impurity scattering matrix element with the 
isotropic screening model becomes [see Eq. ©J 

Vlr^\Jdze'^'^mz)\ 2 - (Al) 

The Taylor expansion of the exponential provides 

Vt P = j -C i + 0(X). (A2) 

where B — — e 2 /2Aeo£ r , and C\ are numerical constants. 
Since B is independent of the state i, the square of the 



matrix elements including the correlation effect becomes 
\V U ~V 3J \ 2 = \C l -C ] \ 2 + 0{\). (A3) 
and without correlation effects 

Vi + V£= 2 -^ + + 0(1). (A4) 

In the Debye approximation, A cx 1/vT, resulting in 
the different temperature dependence, 

|7«-%| 2 cxl + 0(l/VT) 

(Aoj 

Vi + VlcxT + OiVT), 

clearly showing the weaker temperature dependence 
when including correlation effects. 
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